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Abstract 

A brief overview of some recent publications related to the evaluation of two-loop 
Feynman diagrams is given. 



1 Talk given at the International Symposium on Radiative Corrections CRAD96 (Cracow, Poland, 
1-5 August 1996). To be published in the proceedings {Acta Physica Polonica). 



We shall discuss, very briefly, recent progress in evaluating certain types of two-loop 
Feynman diagrams and related issues. Some emphasis will be put on the activities I was in- 
volved in, in collaboration with F.A. Berends, V.A. Smirnov, J.B. Tausk and N.I. Ussyuk- 
ina. Mainly, a collection of "pointers" to the papers describing the relevant methods and 
algorithms is given. The lack of space makes it impossible to include references to all 
papers containing various applications of two-loop calculations, as well as some problems 
related to on-shell calculations with massive particles. 

As a rule, the physical case of four dimensions is understood. In case of need, the 
dimensional regularization M is used as a regulator. 



1 Two-loop self-energy diagrams with masses 

For some special cases, the results were known long time ago. For example, two-loop 
diagrams contributing to the photon polarization operator in QED were first calculated 
in 1955 [0 (see also in [[J). Some other special cases were considered in refs. M. In all 
examples, some of the internal lines were massless and there were just one or two different 
non-zero masses. As a rule, the results were expressible in terms of trilogarithms Li3, 
except for one of the QCD contributions to the quark selfenergy. 

Indeed, the situation becomes more complicated if one is interested in a diagram 
involving a three-particle threshold with all massive particles. In this case, there are 
arguments || that, for a general external momentum k, the result may not be expressible 
in terms of polylogarithms. The simplest example is the so-called "sunset" diagram with 
three massive propagators. Recently, it was considered in a number of papers || 0. 
The results were expressed either in terms of multiple hypergeometric series or via one- 
dimensional integrals. 

For more complicated diagrams, e.g. for the general two-loop self-energy diagram with 
different masses, the representations in terms of known functions are not available^]. One of 
the ways is to construct integral representations and then calculate the result numerically. 
Various approaches to this problem, including some useful integral representations, were 
discussed in refs. [L0|, [U], [12], 0. Tensor reduction of two-loop self-energy diagrams was 
discussed in |13|]. The problem can be reduced to evaluation of the scalar integrals. 



An analytic approach to the calculation of two-loop diagrams with different masses 
consists in constructing the expansions of these diagrams in different regions. For example, 
when there is no threshold at k 2 = 0, the small momentum expansion is basically an 
ordinary Taylor expansion. For two-loop diagrams, the general algorithm for constructing 



the coefficients of such an expansion was presented in |1J]. In general, the expansion 
converges when k 2 is below the first physical threshold. The coefficients of the expansion 
can be presented in terms of two-loop vacuum diagrams. Some results for these diagrams 



can also be found in pT5|] , whereas the results for higher powers of propagators can be 
obtained using the integration-by-parts technique |16| (this procedure is also described 
in [14]]). The general problem of tensor decomposition of two- loop vacuum diagrams was 



discussed in [|T7| [TE|] . In ref . |H| , conformal mapping and Pade approximations were used 

1 In three dimensions, an essential progress has been recently made in ref. ||. However, the three- 
dimensional case is simpler, due to a very simple (exponential) form of the massive propagator in the 
coordinate space. 
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to evaluate numerical values beyond the threshold (s). In ref. pO|1 , a modified scheme for 
calculating the coefficients of the expansion was proposed. 

The problem of constructing the large momentum expansion of two-loop self-energy 
diagrams (when k 2 is larger than all physical thresholds) was considered in [plj| . To do 
this, the general theory of asymptotic expansions of Feynman diagrams |22| (see also in 
[p3| ) was employed. In four dimensions, the coefficients of the large momentum expansion 
of the master two- loop self-energy diagram involve In (A; 2 ) and ln 2 (/c 2 ). The lowest term of 
the expansion of the master diagram is nothing but the massless integral [24] proportional 
to 6C(3). 

Then, the next question is what to do if one is interested in the threshold behaviour. 
Dealing with the thresholds may be unavoidable even in the small momentum expan- 
sion, when the lowest physical threshold vanishes. In this case, the Taylor expansion 
does not work, and one should use a procedure based on p2| . These algorithms for the 
zero-threshold expansion were described in (see also in |p26[) where all possible zero- 
threshold configurations were considered. Moreover, when there is a large mass parameter, 
one can also use the large mass expansion [22| to describe the non-zero threshold behaviour 
at the small (as compared with the large mass scale) non-zero thresholds. For two-loop 
selfenergies, this procedure was considered in ref. |27|. One should not put any condi- 
tions on the relative values of k 2 and small masses. In this case, two-particle-threshold 
irregularities are in fact described by the one-loop two-point diagrams. 



2 Massive and massless two-loop three-point func- 
tions 

The problem of calculating the two-loop three-point functions is, in general, more compli- 
cated than the two-point case. In particular, there are more external invariants {pl,P2 and 
£>!), an d the structure of singularities is more involved. Furthermore, there are two basic 
topologies, the planar one and the non-planar one, and there is a problem of so-called 
irreducible scalar numerators. 

On one hand, for three-point functions one can still use approaches based on numerical 
integration of parametric integrals, see e.g. in [28, 10, 11, 0. The approach used in [1 I 



to calculate the planar vertex function was recently extended to the non-planar case [^9 

On the other hand, it is also possible to extend the analytic approach based on ex- 
pansion of three-point functions in different regions. For example, when all three external 
momenta squared are small (less than the corresponding thresholds), one can expand the 
three-point function in a triple Taylor series, using the algorithms similar to ones from 



14 1 . Such an approach was developed in |1| where also conformal mapping and Pade 



approximations technique were used. An explicit expression for the projectors yielding 
the coefficients of this triple series was given in [ |T8|[ . 

Some special cases of three-point functions involving zero thresholds were considered 
in |50|| . As in the two-point case [E5|, the general theory of asymptotic expansions \T2\ 



was applied. Analogous technique can be also used in the case when all three external 
momenta squared are above the corresponding thresholds. In particular, one needs the 
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results for the corresponding (in general, off-shell) integrals with massless internal linesQ. 
The exact result for the off-shell massless planar three-point function was obtained in 



ref. ||32|| . In the derivation the "uniqueness" technique j33| was used. This approach was 
also generalized to the case of ladder diagrams with an arbitrary number of rungs, see in 
[ |32| For the two-loop diagram, the result was presented in terms of polylogarithms 



up to the fourth order, Li 4 . In ref. the non-planar diagram was also calculated^, as 



well as the cases when some of the propagators are shrunk. The result for the non-planar 
diagram was presented in terms of the square of an expression involving dilogarithms. 
The problem of irreducible numerators (i.e. scalar numerators which cannot be cancelled 
against the denominators) was also considered in [^5| and some results for the integrals 



with such numerators were obtained. A complete solution to this problem should also 
include an efficient algorithm for calculating the integrals with higher integer powers of 
the propagators and irreducible numerators. 
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